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We consider a dilute gas of dipole moments in an arbitrary harmonic trap and treat both the 
short-range, isotropic delta-interaction and the long-range, anisotropic dipole-dipole interaction per- 
turbatively. With this we calculate the leading shift of the critical temperature with respect to that 
of an ideal gas as a function of the relative orientation of the dipole moments with respect to the 
harmonic trap axes. In particular, we determine those magic angles, where the dipolar shift of 
the Bose-Einstein condensation temperature vanishes. Furthermore, we show for the parameters of 
the ongoing ^^Cr-experiment in Stuttgart that this dipolar shift can be enhanced by increasing the 
number of particles, the geometrical mean trap frequency, and the anisotropy of the trap. 
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I. INTRODUCTION 

Ultracold atomic quantum gases are many-body systems where a wide variety of macroscopic quantum phenomena 
is observable [1-3]. For the original Bose-Einstein condensates (BECs) of alkali atoms, it has been sufficient to describe 
the dominant two-particle interaction by a short-range and isotropic contact potential 

v^r(x-x')=^^(x-x'), (1) 

where M stands for the mass of the particles. The s-wave scattering length as is tuned by using a so-called Feshbach 
resonance [4], where it can be varied over a broad range via an external magnetic field. Recently, a BEG has also 
been realized in a dipolar quantum gas of ^^Gr atoms [5] where the magnetic dipole moments are around six times 
larger than those of alkali atoms. Therefore, the physical properties of such a chromium BEG also depend on a long- 
range and anisotropic magnetic dipole-dipole interaction. Other many-body systems with dipolar interactions are, 
for instance, Rydberg atoms [6, 7] or atomic condensates where a strong electric field induces electric dipole moments 
of the order of 10^^ Debye [8]. Permanent dipole moments in heteronuclear molecules are much larger with typical 
values of 1 Debye, so their dipolar effects could be a few hundred times stronger than those of chromium atoms [9, 10]. 
Such a gas of ultracold heteronuclear molecules is produced either by sophisticated cooling and trapping techniques 
[11-14] or by photoassociation [15-18]. For all those systems the specific aniosotropic dipole-dipole interaction reads 

^DD xj- |x-x'|3j+ 3 ^ >' ^ > 

where m denotes the magnetic (electric) dipole moment and fiQ stands for the permeability constant (the reciprocal 
permittivity constant). Note that this dipole-dipole interaction contains, apart from the usual first term, an 
additional contact term, which renormalizes the divergence at the origin. The latter contribution has been overlooked 
in most textbooks of electrodynamics, but is taken into account in more accurate representations as, for instance, in 
Ref. [19]. The dipolar interaction strength can be tuned for induced dipole moments by varying the field strength 
and for permanent dipoles by using rotating magnetic (electric) fields [20]. Gombining this rotation technique with 
Feshbach resonances, will allow in the near future experiments where the interaction varies from purely contact to 
purely dipolar. The anisotropic dipole-dipole interaction gives rise to new condensate properties as, for instance. 
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FIG. 1: Axes of magnetic dipole-dipole interaction (2) and harmonic trap (3) with frequencies u\ > UJ2 > uis- 

a characteristic anisotropic deformation of the expanding BEC which has recently been resolved experimentally in 
Ref. [21]. Also other interesting dipolar phenomena are accessible experimentally as, for instance, the occurrence 
of a Maxon-Roton in the excitation of dipolar BEC [22] or the instability of the ground state of dipolar BECs [6, 23-25]. 

In the present article we investigate how the critical temperature of a BEC is shifted due to the dipole-dipole 
interaction (2). To this end we consider a dilute gas trapped in an arbitrary harmonic potential 



where uji, uj2, denote the respective trap frequencies. Furthermore, we assume that, due to an additional external 
field, the dipole moments m of all constituents are oriented along one axis, i.e. they are uniformly described by 
m = m (sin a cos (/), sin a sin cos a). We expect that the critical temperature depends crucially on the spherical 
angles a and cf) which characterize the orientation of the dipole moments with respect to the harmonic trap axes. The 
temperature shift should have local extrema when the dipole moments are parallel to one of the harmonic trap axes 
as is illustrated in Figure 1. For instance, if the frequencies are ordered according to oji > 0J2 > oJz, this leads to an 
additional interaction which is attractive along the 3-axis, repulsive along the 1-axis and can be both attractive or re- 
pulsive along the 2-axis. Therefore, the particle density in the harmonic trap is increased (decreased) in case 3 (case 1) 
in comparison with a pure contact interaction. Thus, we expect a resulting positive (negative) shift of the critical 
temperature due to the magnetic dipole-dipole interaction. Furthermore, we follow Ref. [26] and suggest to determine 
differences of the critical temperatures in the three configurations, as then the influence of the isotropic contact 
interaction cancels. Consequently, those differences represent a clear signal of the anisotropic dipole-dipole interaction. 

We start our detailed analysis by briefiy reviewing the calculation of the Bose-Einstein condensation temperature 
for a dilute ideal Bose gas in Section II. Subsequently, we treat the influence of both interactions (1) and (2) in 
lowest order perturbatively. This is justified as the confinement of the harmonic trapping potential (3) removes 
critical long-wavelength fluctuations and reduces the fraction of atoms taking part in nonperturbative physics at the 
transition point [27]. Therefore, we determine in Section III within first-order perturbation theory how the grand- 
canonical free energy and, thus, the particle number equation change under the influence of an arbitrary two-particle 
interaction potential. Additionally, we calculate in lowest perturbative order the location of the critical point where 
the phase transition from a Bose gas to a BEC occurs. To this end we extract the critical chemical potential from 
the first-order contribution to the self-energy. This procedure is, in principle, self-consistent and corresponds to the 
Hartree-Fock mean-field theory. However, we determine the critical chemical potential only in the leading order and 
use it to obtain the critical temperture shift. Note that this provides a reliable prediction since it is only linear in 
the interaction strength. Due to the Ginzburg criterion critical fluctuations in the vicinity of the critial point violate 
the mean-field result in higher than first order of the interaction strength [28] . In Section IV we specialize our result 
to the model interaction potential (1) and (2) and combine it in Section V in order to derive the leading shift of the 
critical temperature. Furthermore, we show how it is possible to tune the strength of this dipolar shift. To this end 
we investigate in Section VI the magic angles where the dipolar effects vanish. Finally, we discuss in Section VII our 
theoretical findings for the special parameters of the ongoing chromium experiment in Stuttgart. In particular, we 
estimate which experimental parameters allow to enhance the differences of the critical temperatures for two of the 
three configurations of Figure 1. 
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II. INTERACTION-FREE CASE 



The quantum statistical properties of a dilute ideal Bose gas are determined by the grand-canonical partition 
function which follows from the functional integral 

= ^ vr 2?V' e--^"" I'^* '^^Z" . (4) 

Here f 'Dip*'Dtp sums over all bosonical field configurations ip(x,T) and V'*(x, r) which are periodic in imaginary time 
T. Furthermore, the undisturbed euclidean action in (4) is 

A^°\i}*,ii]=h j dr J dr'y d^^y d3a;'V*(x,r)G(°)-i(x,r;x',r')^(x',T') (5) 
with the integral kernel 

-1(x, r; x', r') = ^ '^(x - x')5(r -r')^^h^-^A + y(x) - ^ | , (6) 

where V{x) is the external trap potential and /U stands for the chemical potential. The grand-canonical free energy of 
such an ideal Bose gas is given by !F^°^ = —{lnZ^^^)//3 and may be evaluated from (4) as J^^^^ = Tr In G^°^~^/p. For 

the harmonic trap potential (3), usually the condition hPu) 1 holds, where a) = {lui uj2 ^s)^^^ denotes the geometric 
mean of the trap frequencies. This allows to determine the grand-canonical free energy within a semiclassical treatment 
[29] 

where (^a{z) = X^^i z^/n"' represents the polylogarithmic function. Now we are interested in thermodynamic prop- 
erties which result from a fixed average particle number A'' = —dJ^^°^ /dfi: 

In particular, we are interested in calculating the critical temperature where a macroscopic occupation of the groimd 
state sets in. The critical point can be immediately read off from Eq. (8) to be fi'^^ = 0, since there the polylogarithmic 
function starts to diverge. However, this critical point, where the phase transition from a Bose gas to a BEC occurs, 
also coincides with a divergence of the correlation function according to the theory of critical phenomena [30]. To this 
end we determine the correlation function of an ideal Bose gas from the function integral 

GW(x,r;x',r') = j V(x, t)V*(x', r') e--^'"'['^*''^l/^ (9) 

The semiclassical result reads 

GW(x,r;x',r') = j e'^^--^'^'^ ^^^{t - T')f^^e-^[^^'' i"^)^ 

n=l 



(10) 



(11) 



The Fourier-Matsubara transform for such a semiclassical expression is given by 

G^°\p,iJm;^) = £\Te'^'-^ J rf^a;'e-'P-'/ftGW (^x + y , r; x - ^, 
where uim = 2TTm/h(3 denotes the Matsubara frequencies. Inserting (10) in (11) yields 

GW(p,a;„;x)= — , (12) 
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which reveals expHcitly that (9) represents the functional inverse of the integral kernel (6). At the critical point, 
where the correlation function (9) diverges, the integral kernel vanishes. Within the semiclassical approximation of 
the Fourier-Matsubara transform of (6) we thus obtain 



G(°)-i(p,a;„;x)= i + ^1 + y(x) - = . 



(13) 



This equation can only be fulfilled for vanishing momentum p = and Matsubara frequency oJm = at the critical 
chemical potential 

=miny(x). (14) 

In case of the harmonic potential (3) the minimum occurs in the center of the trap. Thus, we reproduce the critical 

chemical potential /ic"'' = 0, which corresponds to the semiclassical value of the ground-state energy. Inserting this 
into Eq. (8), we obtain for the transition temperature for an interaction- free Bose gas 

N 



kB 



C(3) 



1/3 

(15) 



where C(a) = X^i^Li 1/^^" is the so-called Riemann zeta function. Note that in more complicated traps, which arise, 
for instance, when the Bose gas is rotated, the potential minimum could occur far away from the origin, thus resulting 
in a non- vanishing critical chemical potential [31]. 

From the critical temperature (15) we read off that hp^^uj ~ N~^/^, i.e. the condition for the semiclassical treatment 
is fulfilled within the thermodynamic limit N oo. For a finite number of particles A'', it is necessary to go beyond 
the semiclassical approximation and to take into account quantum corrections in a systematic way [32-35] . In lowest 
order this finite-size effect in harmonic traps follows from identifying /Xc"'' not with the semiclassical value of the 
ground-state energy but with its quantum mechanical value Eg = 3?EJ/2, where ZJ = (wi + (jJ2 + W3)/3 denotes the 
arithmetic mean of the trap frequencies. Evaluating the polylogarithmic function in (8) perturbatively up to the 
leading order, we obtain for the critical temperature the finite-size effect 



VTi^Vjrs 2C2/3(3)wAfV3 ■ 
The reduction of the critical temperature of an ideal gas decreases with increasing the particle number N. 



III. FIRST-ORDER PERTURBATION THEORY 



In this section we investigate the influence of an arbitrary weak 2-particle interaction upon the thermodynamic 
properties of dilute Bose gases within Feynman's diagrammatic technique of many-body theory [36-41]. 



A. Feynman Rules 

We start with defining the grand-canonical partition function of the full problem 




and the associated correlation function 

G(x,r;x',r') = | j> V(x, r)V*(x', r') e-(-^™[**'*l+-^""'[**'*l)A . (ig) 

In addition to the interaction-free contribution to the action (5) we have also to take into account the interaction 
contribution. Due to the diluteness of the Bose gas we can restrict ourselves to the effects of a two-particle interaction 
y(i„t)(x-x'): 

^('°*)[V'*,^] = \j^ ^ dT J (fx J d^x' y(^"*)(x-x')V*(x,T)V*(x',r)V(x,T)V(x',r). (19) 
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By expanding the functional interals (17), (18) in powers of y('"*^(x — x'), we obtain expressions which consist of 
interaction-free expectation values. These are evaluated with the help of Wick's rule as a sum of Feynman integrals, 
which are pictured as diagrams constructed from lines and vertices. A straight line with an arrow represents the 
interaction- free correlation function (9): 

x',r' = gW(x,t;x',t'). (20) 

Furthermore, spatio-temporal integrals over the two-particle interaction potential are pictured by two vertices con- 
nected by a dashed line 

^"4' = -^J^^dT j (fx y" dV - x') . (21) 

In the following we apply these Feynman rules in order to determine the grand-canonical free energy and the self- 
energy within first-order perturbation theory for the harmonic trap potential (3) and an arbitrary interaction potential 

y(int)(-x_x'). 



B. Grand-Canonical Free Energy 

Up to the first order in the two-particle interaction the grand-canonical free energy = — {\q.Z)II3 reads with the 
diagrammatic rules mentioned above 



(22) 



By passing we note that those two and all higher-order connected vacuum diagrams follow together with their proper 
weights from solving a graphical recursion relation [42]. The first term in (22) is the interaction- free contribution 
to the grand-canonical free energy (7). The second and third term in (22) are called the direct or Hartree-like and 
exchange or Fock-like vacuum diagrams, respectively, which correspond to the following analytical expressions: 

= ^£^dT j (fx y"dVy('"*)(x-x')G(°)(x,T;x,T)G(°)(x',T;x',T), (23) 



= ^ dT j (fx y"ciVF('"*Hx-x')G(°)(x,T;x',r)G(°)(x',r;x,T). 



(24) 



Both contain the interaction-free correlation function (10) with equal imaginary times. In order to guarantee the nor- 
mal operator ordering within the functional integral formalism, this equal-time correlation function G*^°''(x, t;x',t) 
must be interpreted as G^^^ (x, r; x', t+) [38, 39]. Here we have introduced r+ as an imaginary time which is infinites- 
imally later than r. With this and a Fourier transformation of the interaction potential we obtain for the harmonic 
trap potential (3): 

^ 2{h(3Qf i^J^^ {n + n'f J (2^ ^''"'Hq) exp + «')/" ^^^^ 



C. Self-Energy 

In order to determine the location of the phase transition from a Bose gas to a BEC, we follow the physical notion 

elucidated in the ideal case of the previous section and investigate when the correlation fimction (18) diverges, i.e., 
when its functional inverse G~^(x, r; x', r') vanishes. Due to the two-particle interaction, it decomposes into the 
integral kernel (6) and the self-energy S(x, t;x',t') according to 



G-i(x, r; x', r') = G^") -^(x, r; x', r') - S(x, r; x', r') , 



(27) 
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which is equivalent to the Dyson equation [36-41] 
G(x,t;x',t') = GW(x,t;x',t') 

d^'x" j d^x'" j dr" j dT"'G(x,T;x",r")S(x",T";x"',T"')G(°)(x"',r"';x',T'). (28) 

Therefore, a perturbative evahiation of the self-energy yields via (28) self-consistent interaction corrections to the free 
correlation function (9). In lowest perturbative order the self-energy is defined by the diagrams 

S(x,t;x',t') = O + +••••, (29) 

x,r^W x',r' 

which correspond to the analytical expressions 

E(D)(x,r;x',T') = ^ 5(r - r')5(x - x') y" dV W'"*) (x - x")G(°) (x", r; x", r) , (30) 
E('^)(x,T;x',r') = ^ (5(r - tOV^'"*^^ - x')G(°)(x, r; x', r) . (31) 

Note that the diagrams (29) follow from amputating a line in the connected vacuum diagrams of the grand-canonical 

free energy [43, 44]. Thus, using the first-order perturbative approximation (29) for the self-energy corresponds via 
the Dyson equation (28) to a nonperturbative Hartree-Fock mean-field theory for the correlation function. 

To find the zero of (27), we consider its Fourier-Matsubara transform. Together with (6) this leads to the equation 

where the Fourier-Matsubara transform of the self-energy follows in analogy to Eq. (11) from 

E{p,u;rn-,^) = £''dTe''--- J d^'x' e-'P^'/'^^ (k+^,t-,^-^,0^ . (33) 

Eq. (32) is fulfilled up to the lowest perturbative order at vanishing momentum p = and Matsubara frequency 
Wm = as in the previous section. Furthermore, we conclude that the additional contributions (30), (31) in (32) do 
not change the location x^in = of the potential minimum. Therefore, the critical chemical potential is given up to 
first order by = — fiS(O,0;O). Identifying the direct and the exchange contribution with ^^P'' = — fi,S''^-'(0, 0; 0) 
and /x^^^ = — fiE^^^(0, 0; 0), we obtain fic = mI^^ + /^c^^ + • • •, where the respective terms read according to (10) 



IV. MODEL INTERACTION 



In this section we specialize our general formulas for the two-particle interaction between dipolar bosons which 
contains both the contact interaction (1) and the dipole-dipole interaction (2). At first, we briefly discuss whether 
applying this model interaction is physically reasonable. Recently, it has been suggested in Refs. [45, 46] that the 
s-wave scattering length could strongly depend on the dipole moment. However, there it has also been shown 
that for dipolar interaction strengths, which are not larger than the s-wave scattering strength, the latter is only 
rescaled by a moderate factor and remains positive. For the calculations in the present work we assume that this 
condition is fulfilled. Furthermore, we mention that the contact interaction (1) represents an effective pseudopotential, 
see e.g. Refs. [36, 47]. An analogous pseudopotential for anisotropic interactions has recently been introduced in 
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Refs. [48, 49] which contains, apart from the nsual dipolc-dipolc interaction (2) in the s-d scattering channel, also 
another part from the d-s scattering channel. The contribution of the latter turns out to be nonlocal in momentum 
space and is therefore difficult to handle. Here we assume that we are not in the vicinity of any dipolar shape resonance 
and, hence, adopt the anisotropic interaction in the original form (2). Since there exists up to now no experimental 
evidence that a scattering in the d-s channel is relevant, our model seems to be valid for all current experimental 
situations. 



A. Fourier Representation 

Now we combine both interactions (1), (2) and write our model interaction in the following form: 



l2 



M ' ' 47r 1 |x-x'|5 |x-x'|3 
Here we used the effective scattering length 



(36) 



for describing the strength of the contact interaction. Note that both terms in the effective scattering length (37), i.e. 
the original s-wave scattering length as and the magnetic dipole contribution, are physically inseparable. Therefore, 
we have to identify (37) with the experimentally measurable scattering length of atoms with magnetic dipole moments. 
In the following we will use the dimensionless parameter 

which is a measure of the strength of the dipole-dipole interaction relative to the effective scattering energy. In 
the limit of a vanishing scattering length ag ^ the model interaction (36) describes the pure dipolar interaction 
(2). This corresponds to the value cdd = 1, which plays a special role by investigating the stability of a dipolar 
BEC at zero temperature within the Thomas-Fermi approximation [24, 25]. Under certain conditions one expects at 
this dipole interaction strength that a cloud of dipolar particles collapses into a thin wire along the magnetization 
direction. In contrast, the value cdd ^ oo is only achieved in the limit Us — fiQirn^ M / 12it?i?' which does not 
correspond to any distinguished situation. Furthermore, if we consider the limit of a vanishing dipole moment m — > 0, 
the model interaction reduces to the pure contact interaction (1). 

Now we determine the Fourier representation of our model interaction (36). The Fourier transform of the contact 

interaction is simply a constant. The Fourier transformation of the anisotropic term has to be evaluated with special 
care, as the corresponding Fourier integral is UV-divergent. One possibility to regularize this singularity is to introduce 
an UV cut-off distance which physically takes into account that the atoms cannot overlap [9]. Thus, the Fourier 
transform is determined for a finite UV cut-ofi^ distance which is, finally, allowed to go to zero. Using spherical 
coordinates, where the dipole moment is described by m = m (sin a cos sin a sin 0, cos a) , we obtain the following 
Fourier transform: 



sin a (cos (p Qi + sin (j) q^) + cos a 



M 3 ■ q2 

-|- 2sinacosa (cos^giga -|- sin(^(j'293) + 2sin^asin^cos(^(j'ig2 • (39) 

Another approach for obtaining this result (39) has been explored in quantum electrodynamics in the context of 
calculating the transverse delta function with a regularization function [50]. Note that the result (39) also follows 
from applying the distributional identity [25] 

1x15^ =Y^.-^(x) + ^^, (40) 



which can be regarded as defining the dipole-dipole interaction as second partial derivatives of the Coulomb potential. 
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B. Evaluation of Integrals 

Now we evaluate the influence of the interaction (39) upon the thermodynamic quantities of interest. To this end 
we remark that the remaining integrals (25), (26) and (34), (35) are of the form 



/(i. 



J-ooV^ J-ooV^ J-ooV^ V «i «i 4J 



<i2 03 / iiorrfiaia2a3 



In order to calculate this quantity, we define the integrals 



/(i) 



J -00 J -00 J -00 \^ ^^DD V «1 «2 «3/ 



0-2 as / ai 012(13 



which occur for the configurations j = 1, 2, and 3 of Figure 1, respectively. Due to the symmetry of its integrand, 
the full problem (41) is then solved by 



sin2acos2</> /(D ( ^, ^) + sin^asin^^ ( ^, ^) +cos'a f^,^) . (43) 



j(int) / ^ ^ 

Va2' as 

Prom this representation we obtain immediately that 7^'"*)(ai/a2, ai/as) is extremal for arbitrary parameters ai, a2 
and as only if 



j-(int) ^ ^ j(l) ^ ^ 

V 02 ' as / V a2 ' aa 



a = 7r/2 , (/) = 0,7r : 

a = ./2, <A = ./2,3./2: /O-) ( ^, ^) = ( ^, , (44) 

a = 0, 4>- arbitrary : f ^, = f ^ 

Va2 as/ Va2 as 

Whether these extrema lead to local maxima or minima of depends on the concrete values of the coefficients 
/(I) ^ 7(2) a,nd I^^^ . To calculate them we mention the properties 

J(2) (^^^)^ J(l) f ^, ^ V /(3) f ^, ^) = Jd) f^,^), (45) 
Va2 as/ V«3 aij \a2 azj \ai 02/ 

which simply follow from definition (42). Consequently, we only need to determine one of the three integrals, for 
instance, that for j = 1. Applying the Schwinger proper-time representation l/q^ = dTe~^ ^ [30], we arrive at 
the expression 



\a2 as J eoD 2010203 Jo 



2^l/2 ■ 



(46) 



With the help of substitutions the remaining integral is reduced to incomplete elliptic functions of the first and second 
kind [51, (8.111)] 

m k)= dx , E{ct>, k)= dx -L^ . (47) 

Jo [1 - x^) ' {\ — k^x^) ' Jo [1 — x^) ' 

We present the result in the following form; 

/^^)fe^) = /fe^)-— , (48) 

\a2 as/ \a2 as/ edd 

where the anisotropy function / has to be evaluated separately for the three cases ai < 02, 03 and 02 < Oi, 03 as well 
as as < ai, 02. 
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FIG. 2: Logarithmic plot of the anisotropy function f{rj, k) in Eqs. (49) and (50). Black curves represent the special cases (53) 
with (54). The cross line corresponds to fs{n) = /(k, k) and the lines parallel to the axes depict /(I, k) and /(k, 1). 



C. Anisotropy Function 

At first we consider the case ai < a2, as and perform the substitution t = (l/a^ — l/os) /^^ ^ l/'^i which results in 



/(^, .) = 1 + . L f arcsinv/T^, F f arcsinTT^, 

For the case < a2 < ai we rather use the substitution t = (l/af — l/a^) /^^ ^ ^/'^'i^ yield 
^ 1 + V 3At77 / ^ ( 

fiv,-) - - (i-,2)(.2-i) V^^^-'?'^! ™ 



for K, ?7 < 1 . (49) 




^ :F( arcsin^^^,^^^ I S for 1< 7? < k . (50) 



Finally, for the case a2 < 03 < ai, we have a situation which is analogous to the second case when performing the 
replacement r/ ^ k on the right hand side of (50). Note that all these different representations of the function / are 
equivalent to each other by analytical continuation due to the property [51, (8.127)]. The importance of this anisotropy 
function for dipolar Bose gases has recently been recognized in the context of expanding chromium condensates [52]. 
There the anisotropy function appears within the Thomas-Fermi solution of the underlying Gross-Pitaevskii equation 
at zero temperature. Although our investigation of the critical temperature of dipolar gases represents a complete 
different physical situation, the anisotropy function plays here a similar important role. Figure 2 shows how the shape 
of this anisotropy function depends on its parameters r\ and k. For our further investigations we have to summarize 
the properties of /(?], k), see for comparison also Ref. [52]. At first, we mention its symmetry property 

^{^,K) ^ (51) 

and, secondly, the sum rule 

/('?,«)+/f-,-)+/f-,^)=0. (52) 

Furthermore, we need the analytical continuation of this function into the semi-symmetric regions of the trap potential 
(3) 



/(«,/«) =-2/(l,l/«) = /,(«), 



(53) 
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where the anisotropy function 



2k2 



1 



1 - k2 (1 _ ^2)3/2 



2 fi;2 + 1 3k2 



artanh — k? , < k < 1 

, K = l 

arctan \/ n"^ — 1 , k > 1 



(54) 



t l-K? (k2- 1)3/2 

describes a cylindrically symmetric dipolar BEC Rcfs. [24-26, 53]. This function tends asymptotically to —2 and 1 in 
the limit n ^ oo and k ^ 0, respectively (see Figure 2). 



D. Results 



With the abbreviation (42) we find for our perturbative results (25), (26) and (34), (35) in the three extremal 
configurations j = 1,2,3 of (44): 

^(D) = - /O) f 5^1, C3 3 3 {e^>^) , (55) 



(E) ^ ^0^2 

3A3 



^^^^(1,1)C(0 • (58) 



Here A denotes the thermodynamic de Broglie wave length 



and the generalization of the polylogarithmic function of Eq. (7) is defined according to 

Ca,6,c(^) = E E , . (60) 

n=l n—1 ^ ' 

It is worth mentioning that both exchange terms of the grand-canonical energy (56) and the critical chemical potential 
(58) depend only on the contact interaction as the anisotropy function /(1,1) vanishes in (48). Thus, within the 
semiclassical approach, only the direct contributions carrys the anisotropic character of the dipolar interaction. 



V. CRITICAL TEMPERATURE 

The total grand-canonical free energy in the three extremal configurations 1-3 of Figure 1 follows from (22) as the 
sum of (7), (55), and (56): 

C|,|,|(e^^)+.... (61) 
The full problem has to be studied at fixed particle number = —dT/dfi: 

(1,1,1(6^")+..., (62) 

where we have used the identity C| | |(z) — 2Ci,|,|(2;) following from (60). The critical temperature Tc is obtained 

from Eq. (62) in the limit t Mc = A*c°^ + A*c^^ -|- . . . . Recalling (57), (58), we obtain within first-order perturbation 
theory at the critical point 

|,Q,!,|)_,(|)«,},^.., ,3) 



1 



C4 {e^') 



/0)(5^^,5^^]+7(i)(i,i) 



N = 



2,\^{MCof 



/0)(^.^^]+70-)(l,l) 

\UJ2 i^z) 



N 



C(3) 



3A3(fi/3eW)3 



/0)(i^,i^]+/0-)(l,l) 
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where C(a,5, c) = Ca.b.c(l) is a generalization of Ricmanns zeta-function. The first term in Eq. (63) reproduces the 
interaction-free particle number, thus leading to the critical temperature (15). The interaction term changes this, so 
that the resulting shift of the critical temperature in configurations 1-3 reads according to (48), (53), and (54) 



AT, 



/AT, 



(1) 



(2) 



-cs- 



-cg- 



cs_ HoW?M ./wi 
2 127r;i2Af •^U2''^3 
Cs i^om'^M 



(3) 



-cs- 



A, 



(0) 



2 127r/l2A(0) 

Cs fJLQra^M 
2" 127rft2A(°) 



/ 



The dimensionless prefactor cs for the contact interaction has the value 



Cs = 



3C(3) 



C(2)-C 



1 3 3 

2' 2' 2 



3.426 . 



(64) 



(65) 



Now we discuss the physical implications of our first-order perturbative result (64). By setting m = we obtain 
for all configurations 1-3 in Figure 1 the same downwards shift of the critical temperature with the dimensionless 
prefactor (65). This simple result for the isotropic contact interaction was originally derived in Rof. [28] within 
a mean-field approach. Only recently it was confirmed experimentally by investigating the onset of Bose-Einstein 
condensation in the hyperfine ground state of ^^Rb [54] . Furthermore, our result (64) also shows how the shift of the 
critical temperature depends on the anisotropic magnetic dipole-dipole interaction. 



So far we restricted our detailed analysis to those spherical angles (a, cj)) between the symmetry axes of the interaction 
(2) and the trap (3) which coincide with the extremal cases (44). It is straight-forward to extend our result to arbitrary 
angles. For this purpose we use the master integral (41) and its decomposition (43) into the integral (42) for j = 1, 2, 3. 
With this we calculate the following angle-dependent shift of the critical temperature: 



at; 

^(0) 



sin^a 



cos^^ii/ — ,— -hsm^<^/ — ,— -hcos^a/ — ,— >• (66) 



For the special cases ofa = 0, Q! = 7r/2 and 
turn out to be extremal. 



0, 7r/2, TT, 37r/2, Eq. (66) reduces to the previous results (64), which 



VI. MAGIC ANGLES 



Changing the angles (a, (j)) between the symmetry axes of the trap and the magnetization direction allows us to 

tune the dipolar effect between those maximal and minimal values. In particular, there exist so-called magic angles, 
where the dipolar effects vanish. We read off from Eq. (66) that those angles are given for an arbitrary anisotropy 
due to the sum rule (52) by 



».W = »cco. I ± ^si„^, - c»2, /(^. ^)//(2, ^) I . (67) 

Note that the radicand has to be positive, which restricts the range of values for the polar angles 0. In general, 
the possibility to find a magic angle depends crucially on the anisotropy coefficients. Only in the case that 
f{ijJ3/u)i,cos/u2) < — /(t<;i/a;2, wi/ws) is fulfilled, magic angles always exist. 

Now wc present the results for cylindrically symmetric traps. At first, we mention for lui — UJ2 that the magic 
angle (67) reduces to a value, which is independent of the second anisotropy parameter uii/cos. This value is also 
(^^-independent according to (67) with (53) and amounts to ao = arccot(±l/-\/2) = 54.7° or 125.3° [20, 26]. For 
another symmetric case with loi = we find with (53) that /(cja/tJi, Cij3/ti;2) = /(li;i/ci;2, wi/w-j), thus solutions 
of (67) only exist for sin^ ^ > 1/3. The locations of magic spherical angles for both symmetric cases loi = L02 and 
wi = ws are depicted in Figure 3 a). For the third symmetric case u}2 = 1^3 the situation is similar to wi = ^3 as the 
polar angle (p is only shifted by 7r/2. 
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a) 0° 90° 180° 270° 360° b) 0° 90° 180° 270° 360° 



FIG. 3: The location of magic angles (a, (f)), where the dipolar effect vanishes, a) The solid line stands for a cylinder symmetric 
problem with wi = 0)3 and the dashed line for that with symmetry uji = lo^. b) The solid line stands for the anisotropic 
parameters of the Stuttgart experiment and the dashed line for a cylinder symmetric problem with u)\ = UJ2- 



Furthermore, we discuss the situation in a complete anisotropic harmonic trap by the example of the trap frequencies 
of the Stuttgart experiment [5] which are given by wi = 27r • 581 Hz, uj2 = 271 ■ 406 Hz, and = 27r • 138 Hz. Thus, we 
find from (50) the values /(ti;i/u;2, wi/wa) = —0.6252 and f{uJz/u>i,wz/uJ2) = +0.7356. Specializing the general result 
(67) for the Stuttgart experiment leads to 



ao{4>) = arccot|±y 0.150 sin^c/) + 0.850 cos20 j . (68) 

This behaviour is depicted in Figure 3 b), where we see the location of the magic angles compared with the 4>- 
independent value of a problem with a cylinder symmetry around the z-axis. 



VII. STUTTGART EXPERIMENT 



Here we discuss the consequences of our results for the concrete example of the parameters from the ongoing 
experiments on the Bose-Einstein condensation of ^■^Cr-atoms at the University of Stuttgart [5]. There the geometric 
mean frequency is tD = 27r • 319 Hz. The total atom number is = 100 000, yielding an interaction-free critical 
temperature (15) of about Tc'^^ — 670 nK. The finite-size correction of ri*^"* in the Stuttgart experiment follows from 
(16) and amounts to —1.8 %. This can now be compared with a shift of the critical temperature due to the contact 
and the magnetic dipole-dipole interaction. 

The scattering length of the ^^Cr-atoms is given by a ~ 105 with the Bohr radius as [55], so we obtain from 

(38) eoD = 0.144. The thermodynamic de Broglie wave length has the value A^^'' ~ 5 598 as for the parameter 
of Stuttgart experiment. From Eq. (64) we obtain from the contact interaction a downwards shift of the critical 
temperature by —6.4 %. This is now modified by the magnetic dipole-dipole interaction. At first, we consider case 
1 from Figure 1, in which the magnetization is parallel to the trap axes with the largest moment of inertia. There 
we find from (50) the result f{uji/ijj2, ^1/(^3) — —0.6252, which leads via Eq. (64) to an extra downward shift of 
the critical temperature by —0.29 % due to the magnetic dipole-dipole interaction. For the case 2 of Figure 1 we 
find f{LU2/uj3, 1^2/(^1) — —0.1104, which amounts to an additional downward shift of the critical temperature due to 
the magnetic dipole-dipole interaction by —0.05 %. Finally, for case 3 of Figure 1 with the magnetization direction 
parallel to the trap axes with the smallest moment of inertia, we find f{uj3/uji,uj3/uj2) = -1-0.7356 which results in 
an upward shift of the critical temperature due to the dipole-dipole interaction of about -1-0.34 %. Figure 4a) shows 
the resulting total shift of the critical temperature ATc for the ^^Cr gas with respect to the interaction- free critical 
temperature Ti"^ versus the particle number N . Both the finite-size corrections and the contact interaction lead to a 
huge shift, on top of which the small dipolar effect is seen. 

Furthermore, we we follow Ref. [26] and suggest to measure differences between the critical temperatures in the 
three extremal cases of Figure 1. The large shift caused by the finite-size corrections (16) and the contact interaction 
would cancel due to their isotropic character. Thus, the difference between the shifts is exclusively caused by the 
magnetic dipole-dipole interaction. For the total atom number N = 100 000 and the mean frequency uj = 2tt ■ 319 Hz, 
the temperature difference between the case 1 and 3 amounts to a net effect of 0.61 % of the interaction-free critical 
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FIG. 4: a) Shift of the critical temperature ATc with respect to the interaction-free critical temperature Ti"' for a ^'^Cr gas 
in a harmonic trap with frequencies coi — 2-k ■ 581 Hz, u!2 = 2n ■ 406 Hz, and ojs = 27r ■ 138 Hz calculated from (64) and (16) 
without (straight line) and with (dashed lines) magnetic dipole-dipole interaction for the cases 1-3 of Figure 1. b) Differences 
of the critical temperature shifts versus effective parameter uiN^^^ . The respective dots indicate the present parameters of the 
Stuttgart ^■^Cr experiment. 



In(c^iM) ^ -2-5, 






FIG. 5: Shifts of critical temperatures a) ATi^'-ATj^', b) AT^^^-AT^^^ and c) AT^'^'>-AtP with respect to the interaction-free 
critical temperature ri"' for A*' = 10^ ^'^Cr-atoms in a trap with mean frequency ti; = 27r • 319 Hz versus anisotropy parameters 
uji/uj2 and uji/ujs according to (64). The black curves represent cylindrically symmetric cases and the black squares the values 
of the Stuttgart experiment. 



temperature Tc . This dipolar effect increases with the geometric mean frequency uj and the total atom number N 
as [ujN^^^]^^'^ , as is shown in Figure 4b). Another possibility to enhance the difference of the critical temperatures is 
provided by varying the anisotropy strengths uji/ll!2 and uji/u^ of the harmonic trap potential as seen in Figure 5. For 
instance, in order to increase the difference of the critical temperature between the first and third configuration, we 
need to increase the aspect ratio lui/u!2- Furthermore, it is preferable to work with moderate anisotropy parameters 
ujxjio-i < 1 (pancake like) than in the regime uji/lus > 1 (cigar shaped) of the Stuttgart experiment, as is indicated 
by Figure 5 b). 

Finally, we elucidate how to achieve measurable dipolar effects at the critical point. At first, varying the anisotropy 
parameters only allows to increase the critical temperature difference ri^^— Ti^'' with respect to the interaction- free 
value from currently 0.61 % up to a maximal value of around 1.39 %. Note that the latter represents an upper limit 
according to our perturbative calculations. Not so strictly limited is the possibility to amplify the magnetic dipolar 
effect by increasing the particle number or the mean trap frequency. For instance, a doubling of the latter yields by 
constant particle number an increase of the effect by the factor 1.41. An even higher increase of the particle number 
seems to be promising, but is in fact physically limited by three-body losses. 



VIII. CONCLUSION 



So far, the dipolar nature of nowdays available BECs has only been resolved in expansion experiments [21]. The 
analysis of the present article shows that it could become possible to detect a signal of the underlying dipole-dipole 
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interaction also from measuring the critical temperature of the onset of Bosc-Einstein condensation. Most sensitive 
is the difference of the critical temperatures where the magnetization direction is parallel to the axes of the harmonic 
trap with largest and smallest moment of inertia. We have shown quantitatively that this temperature difference 
increases with the number of chromium atoms, the geometrical mean frequency, and the anisotropy of the trap. 
Furthermore, our results will, certainly, be useful also for other dipolar systems with a tunable dipole moment, like 
heteronuclear molecules in low vibrational states [11-13, 15-18], where the dipolar effect will be larger. 
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